Linear in-plane magnetoconductance and spin susceptibility of a 2D electron gas on a 

vicinal silicon surface 
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In this work we have studied the parallel magnetoresistance of a 2DEG near a vicinal silicon 
surface. An unusual, linear magnetoconductance is observed in the fields up to B — 15 T, which 
we explain by the effect of spin polarization on impurity scattering. This linear magnetoresistance 
shows strong anomalies near the boundaries of the minigap in the electron spectrum of the vicinal 
system. 



Over the last few years a number of reports have ap- 
peared on observations of anomalous positive magnetore- 
sistance (negative magnetoconductance, MC) in high- 
mobility 2D electron and hole gases in the field paral- 
lel to the 2D plane-. Such interest is fuelled by the fact 
that this negative MC is directly related to other unusual 
properties of high- mobility systems, such as 'metallic' be- 
haviour and the transition from 'metal' to 'insulator'—. 
However, until now all studies of in-plane MC were de- 
voted to the range of low electron densities, N s , near 
the 'metal-insulator' transition. In the case of a 2D elec- 
tron gas (2DEG) in Si-MOSFET structures the densities 
studied were below 10 12 cm -2 . 

The focus of our study is the in-plane MC in a vicinal 
Si-MOSFET at high electron densities, N s > 10 12 cm" 2 . 
It is well known 2 *&ii& that this system has a superlattice 
potential on the Si surface, which results in a minigap in 
the energy spectrum. Previously, the study of the slow 
electron diffraction by atomically pure vicinal Si surfaces 
cut at small angles 9 to (100) plane (as shown in Fig. la) 
revealed ordered steps, which do not disappear even when 
specimens are heated up to 1100°C, even in the presence 
of hydrogen or oxygen 6 . The size of these steps was in 
agreement with theoretical predictions for perfect high- 
index surfaces^. This periodic structure is considered 
to be the reason for the appearance of the superlattice 
potential^. 

The minigap in the energy spectrum has an inter- valley 
character 4 , which implies a strong inter- valley interac- 
tion when the electron energy approaches the minigap. 
This also leads to the appearance of a logarithmic diver- 
gence in the density of states D(E) at energies close to 
the lower edge of the minigap, Fig. 1(b). D(E) differs 
from zero within the minigap because the superlattice is 
one-dimensional and does not produce a full gap in the 
spectrum. The discontinuity in D(E) also appears at 
the upper edge of the minigap. When the smearing T is 
less than the width of the minigap A, transport coeffi- 
cients exhibit singularities as the Fermi energy Ep passes 
through the boundaries of the minigapi^i 

Such systems with a minigap are of interest because 
they can realise two different situations: with an isotropic 
Fermi surface (FS), and an anisotropic one, dependent on 
the position of the Fermi level. These two cases corre- 
spond to different ranges of electron densities, where the 



properties of the 2DEG significantly differ: N s < N A 
(isotropic FS) and N s > N A (anisotropic FS), where 
N A — 7r(0.15/L) 2 is the density corresponding to the on- 
set of the superlattice minigap. (Here L = a/(2sin(0)) is 
the superlattice period, 9 is the angle between the vicinal 
and (100) silicon surface, a = 5.43 A is the lattice con- 
stant of Si.) The properties of the 2DEG in the first range 
are identical to those of an ordinary (100) Si-MOSFET, 
while in the second range the properties of the 2DEG 
become strongly modified^. 

We experimentally investigate the MC in both ranges 
of N s . A linear magnetoconductance is observed in mag- 
netic field parallel to the plane of the 2DEG, in a surpris- 
ingly large range of fields up to 15 T. To our knowledge, 
this is the first observation of such an effect. We de- 
scribe the linear decrease of the conductance with B\\ 
in terms of the screening model^, where impurity scat- 
tering gets changed when the 2DEG is spin-polarised by 
parallel magnetic field. The performed analysis suggests 
that the slope of the linear MC in the first range of N s 
(isotropic Fermi surface) is in agreement with theoretical 
expectation for the scattering dominated by impurity po- 
tential. In the second range of N s a significant decrease 
of the magnitude of the MC is observed at the bound- 
aries of the minigap. This can be attributed to either the 
modification of the transport properties of the electrons 
near the points of topological transitions of the Fermi 
surface 7,11 , or suppression of the spin susceptibility. 

The studied samples are Hall-bar Si-MOSFETs fabri- 
cated on a vicinal Si(17,2,2) surface, tilted from (100) by 
an angle 9 — 9° 40' around the direction [011]. The su- 
perperiod in this case is L = 16.2 A, which corresponds to 
N A ~ 2.6 x 10 12 cm -2 . As this surface is just slightly dif- 
ferent from the (811) surface at 6 = 10°, one can expect 
qualitatively the same dispersion relation and other prop- 
erties as those for Si(nll) described in detail in Ref.4. We 
study the behaviour of the 2DEG for both orientations of 
the Hall bar - along and perpendicular to the superlat- 
tice axis. In the latter case, shown in Fig. 1(a), the peak 
mobility of the 2DEG is ~ 25000 cm 2 /Vs. Our measure- 
ments have been carried out in magnetic fields up to 15 
T at T ~ 50 mK. The experiments have been performed 
at large electron densities: N s from 1 x 10 12 to 5 x 10 12 



cm 



In magnetic field parallel to the plane of the 2DEG 



2 




3 6 9 12 15 3 6 9 12 15 



B„(T) B M (T) 

FIG. 1: (a) A schematic diagram showing the geometry of the 
studied sample in the case of the Hall-bar oriented perpendic- 
ular to the super-lattice axis, (b) Density of states of 2DEG 
in a vicinal sample with super-lattice potential at the surface. 
Do is unperturbed density of states (as in (100)Si MOSFETs), 
and A is the minigap width, (c-d) Longitudinal conductivity 
versus in-plane magnetic field for different electron densities 
(symbols), with linear fits (solid lines). 



tures is seen in Fig. 2 (a,b) - these features seem to be 
significantly weakened when a strong in-plane magnetic 
field is applied. 

To analyse the linear MC let us first consider the range 
N s < Ng, where, as mentioned above, the properties of 
a 2DEG near the vicinal surface and the (100) surface 
are expected to be equivalent. Also, it is well estab- 
lished that scattering in the low density regime is dom- 
inated by impurity scattering rather than by interface 
roughness^. In this region of N s the slope of the linear 
MC decreases monotonically with increasing density as 
it has been shown above. 

Among several models of negative parallel-field MGi£, 
there are two which indeed predict linear MC: one is 
the interaction theory by Zala, Narozhny and Aleiner— , 
and the other one is the screening model by Dolgopolov 
and Gold£. We start our consideration from the interac- 
tion theory, which addresses the effect of in-plane mag- 
netic field on quantum correction to conductivity due to 
electron-electron (e — e) interactions in 2D disordered sys- 
tems. This theory predicts linear MC in relatively strong 
magnetic fields and in the ballistic regime, when the pa- 
rameter huTr/h ^> 1. According to Ref. [10] the strong 
field criterium (g* fj,BB\\/2ksT 3> 1, where g* - Lande g- 
factor, fiB - the Bohr magneton) is realised in our exper- 
iment already at Bit > 0.3 T. However the values of the 
parameter ksTT/Ti range from 0.003 to 0.015, implying 



we have observed negative MC in the entire range of N s . 
An example is shown in Fig. l(c, d) for the Hall bar 
oriented parallel to the superlattice axis. In the whole 
range of fields the conductivity decreases linearly with 
Bi | . The slope of this linear dependence decreases mono- 
tonically with increasing density up to N s ~ 2.5 x 10 12 
cm -2 , which is close to iV A , Fig. 1(c). However, at larger 
densities the slope starts changing nonmonotonically, as 
seen in Fig. 1(d). It is important to emphasise that MC 
remains linear in the entire density range. 

In Fig. 2 we show the conductivity measured as a 
function of electron density at different magnetic fields 
from to 15 T, for both Hall-bar orientations. In zero 
magnetic field the so called "f2"- and "W"- features are 
clearly seen in a xx (N s ) for the two orientations, respec- 
tively (the names reflect the shapes of the two depen- 
dences in the minigap). These well known features were 
previously observed in vicinal systems^*. They originate 
from the superlattice structure of the vicinal surface and 
indicate the very good quality of our samples. 

One can notice that the effect of magnetic field on the 
conductivity is minimal in the characteristic points cor- 
responding to the lower and upper boundaries of the su- 
perlattice minigap, N A and N A , marked by the vertical 
dotted lines in Fig. 2 (a,b). This is seen for both Hall- 
bar orientations and corresponds to the non-monotonic 
dependence of the linear MC on electron density in Fig. 
1(b). Also, a transformation of the "W" and "f2" fea- 
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FIG. 2: Density dependence of the conductivity at different 
in-plane magnetic fields: B\\ =0, 2, 4, 6, 8, 10, 12, 15 T; (a) 
for the Hall-bar oriented perpendicular to the supperlattice 
axis, (b) along the supperlattice axis. The dotted lines mark 
the boundaries of the superlattice minigap at iV A and iV A 
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that the system is in the diffusive regime, ksTT/h <C 1. 
The estimation of MC caused by electron interactions in 
this regime^ shows that effect of e — e interactions is neg- 
ligible, as the magnitude of MC (proportional to \n(B)) 
raises to only about 1% at Bu ~ 10 T. 

In the screening model^, developed for the system 
with isotropic Fermi surface and at T = 0, the MC 
effect originates from the change in the screening of 
scattering potential caused by the difference between 
Fermi momenta k F - and k F+ for the two spin-split sub- 
bands. It predicts the change of the Drude conductivity 
with the following negative MC in the low-density case, 
q s ^> 2k Fl where q s is the screening wave number. At 
small magnetic fields Bu such that Bu < 0.2Bs, where 
Bs = 2Ep/ g* fiB is the field of the full spin polarisation, 
this MC is linear in the case of short-range impurity scat- 
tering: 
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where the coefficient a varies in the range from 0.87 to 
0.68 for the electron densities from N s — 1 x 10 12 to 
5 x 10 12 cm~ 2 , Ref. [8-9]. As the latter is the range 
of N s of our experiment, we take the average value of 
a = 0.78 for the approximate analysis. 

In our case q s /2k F ~ 2 — 4 which justifies the approx- 
imation of the low density. To justify the applicability of 
the theory developed for T = 0, we only analyse quanti- 
tatively the data at the lowest temperature T = 50 mK. 
Simple estimation using standard parameters for (100) Si 
MOSFETs (g* ~ 2.5, m* = 0.19m e ) gives then 10 - 15% 
magnetoconductivity in magnetic field B ~ 15 T, which 
is close to what is seen in Fig. 1(a). The screening theory 
analyses the effect of magnetic field in both cases of im- 
purity and roughness scattering provided Bu = 3$, Ref. 
[9]. However, the expression for the small field MC has 
been obtained only for the impurity scattering, Eq. (Q). 

As a first approximation, we analyse the slope of the 
linear dependence in the whole range of electron densities 
using Eq. and neglecting the influence of interface- 
roughness scattering. It follows from the above definition 
of Bs that Bg 1 cx g*m* and hence it is proportional to 
the spin susceptibility of the electron gas \- (Indeed, for 
a 2DEG in a Si MOSFET E F = Trh 2 N s /2m* and x = 
2/i 2 B g*m* /ttH 23 ) . The electron density is known from the 
measurements of Shubnikov-de Haas (SdH) oscillations. 
From the slope of the linear MC the product g*m* is 
then found and normalised by <7o m f>, where mt = 0.19m e 
is bulk electron mass (m e is the free electron mass) and 
go = 2.0 is the Lande g-factor in bulk silicon. The data 
for the both Hall-bar orientations are plotted in Fig. 3 (a, 
b). In Fig. 3(a) we also show the g*m* results obtained 
by three different experimental groups using the analysis 
of SdH oscillations in ordinary (100)Si MOSFETa 13 i 14 i 15 . 
At lower densities our results show a monotonic increase 
of spin susceptibility with decreasing N s . This trend at 
N s < 2.3 x 10 12 cm~ 2 as well as the magnitude of g*m* 
agree well with the previous results. The above range is 
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FIG. 3: (a, b) Density dependence of the product g*m* (pro- 
portional to spin susceptibility x) normalised by the product 
of the bulk parameters gorrib: (a) for the Hall bar oriented 
perpendicularly to the superlattice axis, (b) along this axis. 
Data for (100)Si MOSFETs: (o) - from Ref.[13], (•) - from 
Ref. [14], (+) - from Ref. [15]. (c) Density dependence of the 
conductivity at B — for the Hall-bar in case (a) (the same 
data as in Fig. 2 (a). 



exactly the first range of N s we referred to before (N„ < 
N A ), where the Fermi surface of the vicinal sample is 
isotropic, as it is in the case of (100)Si, and where the 
screening theory is valid. This result also indicates that 
in this range of _/V s the effect of the interface-roughness 
scattering is negligible. 

A strong deviation from the monotonic behaviour of 
spin susceptibility is observed at larger densities (N s > 
N A ) if the analysis is carried out using the same ap- 
proach as above. A drop of g*m* by more than a factor 
of two is seen at two different N s for the Hall-bar oriented 
perpendicularly to the superlattice axis, Fig. 3(a). The 
two minima also appear for the other Hall-bar orientation 
(along the superlattice axis), Fig. 3(b), although they are 
less pronounced, presumably because of the smaller mo- 
bility of this sample (the difference in the conductivities 
is seen in Fig. 2). It is clearly seen from the compar- 
ison to Fig. 3(c) that these minima coincide with the 
dips of the "ST and "W" features in a(Ns) - the points 
where the Fermi level crosses the boundaries of the su- 
perlattice minigap. (The pronounced decrease of MC at 
these points is also seen in Fig. 2 a, b.) It is important 
to notice that the changes of the conductivity with N s 
at the minigap boundaries do not exceed 10%, Fig. 2, 
while a much more drastic decrease is seen in the density 
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FIG. 4: (a) Typical dispersion relation for a 2DEG in a Si 
vicinal system—, (b-e) Fermi surfaces at different positions of 
the Fermi level Ef marked by horizontal dotted lines in (a). 



dependence of the spin susceptibility. 

We have to note, however, that the screening theory is 
not expected to be valid in the range of N s > N A . At 
the same time, it is seen in Fig. 3(a) that in the middle 
of the minigap the obtained values of g*m* are close to 
those obtained for (100)Si using SdH oscillations. This 
points out that the anomalies in the MC arise only at the 
boundaries of the minigap. 

Let us discuss possible reasons for the dramatic de- 
crease of the MC at the minigap boundaries. A typical 
dispersion relation of the vicinal system in Si is shown in 
Fig. 4(a). Different configurations of the Fermi surface 



are shown in Fig. 4(b-e) for four positions of the Fermi 
level. It is seen on the plot that the boundaries of the 
minigap in the spectrum correspond to the topological 
transitions of the Fermi surface: firstly the two isotropic 
and independent FSs coalesce (Fig. 4c), and then they 
form two inclosed surfaces (Fig. 4e) . It is known that the 
conductivity and thermopower of the electron system is 
strongly modified at the points of the topological tran- 
sitions (N s — N A , N A )^. In our case we see a strong 
decrease of the linear MC at these points. 

The exact theory of MC for complicated FS is not 
developed. However, if one assumes that the screening 
theory^ can still be applied at the boundaries of the 
minigap, our results indicate that the spin susceptibility 
sharply decreases in the points of the topological tran- 
sitions. Continuing this logic, it is important to note 
that previous investigations^ have shown that the den- 
sity of states at the minigap boundaries changes only 
within 10%, which would correspond to a small change 
in the effective mass. This implies that it is mainly the 
g-factor that is affected by the topological transitions in 
the 2DEG spectrum of the vicinal system. 

An increase of the g-factor, caused by enhanced ex- 
change interaction, is well knowni&. In contrast, in our 
case a significant decrease of g* (by a factor of two) 
is observed. This effect could possibly be caused by 
the strong inter-valley interaction in the vicinal struc- 
ture (giving rise also to the "fi" and "W" features in 
the conductivity^). We have to say, however, that the 
above speculations are based on the assumption that the 
screening theory^ is applicable within the minigap. An 
appropriate theory for such a case does not exist at the 
moment, and we hope that our experiments will stimu- 
late its development. 
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